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Abstract 



We present a high statistics analysis of the pure gauge compact U(l) lattice 
theory using the the world— sheet or Lagrangian loop representation. We have 
applied a simulation method that deals directly with (gauge invariant) integer 
variables on plaquettes. As a result we get a significant amelioration of the 
simulation that allows to work with large lattices avoiding the met aest ability 
problems that appear using the standard Wilson formulation. 



I. INTRODUCTION 



D 
> 
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It is well known that pure U(l) lattice gauge theory exhibits a phase transition separating 
a region where photons remain confined (strong coupling regime) and a Coulombian region 
(weak coupling). Nevertheless, the order of such transition is rather controversial |1]|§. 
Simulations on hypercubical lattices using the Wilson form of the action shown the existence 
of a latent heat that seems to persist in the infinite volume limit |M].If this transition is 
actually of first order, hence the gauge U(l) lattice theory would be lacking of a continuum 
limit. 

The seriousness of such a situation has motivated the appearance of two lines of study. 
The first one is intended to check if the unexpected first order nature of this phase transition 
is an "accident" of the simulation, originated by the imposition of an hypercubical lattice || , 
the periodic boundary conditions ||, etc. The second line of work goes directly to the the 
physical origin of the phase transition, the monopole condensation , and to the analysis of 
all related effects, as the presence of monopole percolation ||, that can shadow the behavior 
of the phase transition. 
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Although the formulation of the compact pure gauge U(l) theory on a lattice is very 
simple, i.e. there are only phases associated to the links, obtaining reliable numerical re- 
sults from simulations is a rather difficult task. The earlier numerical simulations of such 
a theory indicated actually a continuous phase transition PJlQfl. The evidences which fa- 



vor a first order transition appeared only when larger lattices and higher statistics were 
accessible |TT|--[T3| . However, numerical simulations on larger lattices are very difficult. This 
is originated mainly by the lack of tunneling between phases that become incredibly stable 
and difficult the statistical analysis. 

The ordinary Wilson formulation of QED has the links as the basic variables represent- 
ing the gauge field and therefore involves a gauge redundancy. A different implementation 
of QED on the lattice which avoid the gauge redundancy can be attained using a descrip- 
tion directly in terms of the world-sheets of gauge invariant quantities or loop excitations. 
This is the Lagrangian counterpart of the former introduced Hamiltonian loop representa- 
tion [|TJ,|T^]. Although the non-canonical nature of the loop algebra had made elusive this 
Lagrangian formulation, in a previous set of papers [[U||17]] a procedure to cast loops into 
the so-called world-sheet formulation was introduced. 

In this paper we present a large statistics analysis of the pure gauge U(l) lattice the- 
ory using the world-sheet formulation. The lattice action used is equivalent to the Villain 
form. In this formulation basic variables are integer variables on the plaquettes. As we will 
comment, the thermodynamical behavior over the phase transition is different to the one 
obtained from a standard simulation. The persistent metaestability problems that appear 
in the simulations using the traditional formulation are absent in the case of the world-sheet 
formulation . This therefore makes this formulation a useful computational tool. 

The plan of the paper is as follows. In Sec. II we summarize the basics of the world-sheet 
approach. In Sec. Ill we present the simulation method and the results using the world- 
sheet action corresponding to the Villain U(l) theory. Sec. IV is devoted to the analysis of 
the critical exponents of the transition. Finally, sec. V contains the summary of the work 
and the main conclusions. 



II. THE LAGRANGIAN LOOP APPROACH OR THE WORLD-SHEET 

FORMULATION 



QCD is generally defined in terms of local fields, quarks and gluons, but the physical 
excitations are extended composites: mesons and baryons. There is an alternative quantum 
formulation of Yang-Mills theories directly in terms of these extended excitations namely, 



the loop representation [|n],|IJ|. The basis of the loop description of gauge theories can 
be traced to the idea of describing gauge theories explicitly in terms of Wilson loops or 
holonomies [18-20 since Yang |^T| noticed their important role for a complete description 
of gauge theories. Loops replace the information furnished by the vector potential. A 
description in terms of loops or strings, besides the general advantage of only involving the 
gauge invariant physical excitations, is appealing because all the gauge invariant operators 
have a simple geometrical meaning when realized in the loop space. 

The loop based approach of ref. flH| , |To| describes the quantum electrodynamics in terms 
of the gauge invariant holonomy (Wilson loop) 
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W( 1 )=exp[iei A a (y)dy% (1) 

and the conjugate electric field E a (x) . They obey the commutation relations 

[E a (x), #(7)] = e / 5{x - y)dy a W( 1 ). (2) 



These operators act on a state space of Abelian loops ^(7) that may be expressed in 
terms of the transform 

V(7) = / d„[A] < 7 I A >< A |V >= / d M [A]V[A]exp[-ie ^ A a dy ]. (3) 

This loop representation has many appealing features: first, it allows to do away with 
the first class constraints of gauge theories. That is, the Gauss law is automatically satisfied. 
Second, the formalism only involves gauge invariant objects. Third, all the gauge invariant 
operators have a transparent geometrical meaning when they are realized in the loop space. 

When this loop representation is implemented in the lattice it offers a gauge invariant 
description of physical states in terms of kets | C >= W(C) | >, where C labels a closed 
path in the spatial lattice. Eq.(^) becomes 

[E l ,W(C)}=N l (C)W(C), (4) 

where / denotes the links of the lattice, E(l) the lattice electric field operator, W(C) = 
Yli e c U(l) and Ni(C) is the number of times that the link I appears in the closed path C. 
In this loop representation, the Wilson loop acts as the loop creation operator: 

W(C) I C >=| C'-C > . (5) 

The physical meaning of a loop may be deduced from (|) and fl5|), in fact 

Ei\C>= N t (C) I C >, (6) 

which implies that | C > is an eigenstate of the electric field. The corresponding eigenvalue 
is different from zero if the link I belongs to C. Thus C represents a line of electric flux. 
In order to cast the loop representation in Lagrangian form it is convenient to use the 



language of differential forms on the lattice of ref. p2 |. Besides the great simplifications to 
which this formalism leads its advantages consists in the general character of the expressions 
obtained. That is, most of the transformations are independent of the space-time dimension 
or of the rank of the fields. So let us summarize the main concepts and some useful results 
of the differential forms formalism on the lattice. 

A k-form is a function defined on the k-cells of the lattice (k=0 sites, k=l links, k=2 
plaquettes, etc.) over an Abelian group which shall be R, Z, or U(l)=reals module 2ir. 
Integer forms can be considered geometrical objects on the lattice. For instance, a 1-form 
represents a path and the integer value on a link is the number of times that the path 
traverses this link. V is the co-border operator which maps k-forms onto (k+l)-forms. It 
is the gradient operator when acting on scalar functions (0-forms) and it is the rotational 
on vector functions (1-forms). We shall consider the scalar product of p-forms defined 
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< a | (3 >= J2c k a ( c )P( c ) where the sum runs over the k-cells of the lattice. Under this 
product the V operator is adjoint to the border operator d which maps k-forms onto (k- In- 
forms and which corresponds to minus times the usual divergence operator. That is, 

< a | V/3 >=< da\P>, (7) 
<Va\f3 >=< a | d(3 > . (8) 

The co-border V and border d operators verify 

V 2 = 0, d 2 = 0. (9) 
The Laplace-Beltrami operator operator is defined by 

□ = V<9 + <9V. (10) 

It is a symmetric linear operator which commutes with V and d, and differs only by a minus 
sign from the current Laplacian A M A M . From Eq. (|i~0|) is easy to show the Hodge-identity: 

1 = dcr 1 V + vn^a (11) 

After this parenthesis on differential forms on the lattice let us consider the generating 
functional for the Wilson lattice U(l) action: 



rWexp(-f Ecosflp), (12) 

J -IT A p 



A\' 



where the subscripts I and p stand for the lattice links and plaquettes respectively. 
Fourier expanding the expfcos^] we get 



AY 



rWllE^(^)e mA , (13) 

J — 7T 



p n p 

which can be written, using the differential forms language as 

Zw=T, r (^)exp(Eln^(/5))e J <^>. (14) 

In the above expression, Q\ is a real periodic 1-form, that is, a real number 9 G [— ir, tt] defined 
on each link of the lattice; V is the co-border operator; n p are integer 2-forms, defined at 
the lattice plaquettes. By eq. (H) and integrating over 9i we obtain a S(dn p ). Then, the 
partition function can be written as 

Z w cx exp(^ In / np (/?)), (15) 

{n p ;dn p =0} V 

the constraint dn p = means that the sum is restricted to closed 2-forms. Thus, the sum 
runs over collections of plaquettes constituting closed surfaces. 

An alternative and more easy to handle lattice action than the Wilson action is the 
Villain action. The partition function for this action is given by 
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Z v = |(^)^exp(-^||W-27r S || 2 ), (16) 
where 1 1 . . . 1 1 2 =< ...,...>. If we use the Poisson summation formula 

/oo 
-oo 



,g27ri0n 



and we integrate the continuum variables we get 

Z v = (2ttPv)- Np/2 I (^)E ex P(-^- <n,n>+i<n,V9 >), (17) 

where N p in the number of plaquettes of the lattice. Again, we can use the equality: < 
n, V6* >=< dn, 9 > and integrating over 9 we obtain a 5(dn). Hence we get: 

Z v = (2tt/V)-^ /2 E exp(--j- <n,n>), (18) 

{n;dn=0} 

where n are integer 2-forms. Eq. ([Tj^) is obtained from Eq.(|T5|) in the j3 — > oo limit. 

If we consider the intersection of one of such surfaces with a t — constant plane we get 
a loop C t . It is easy to show that the creation operator for this loop is just the creation 
operator for the loop representation i.e. the Wilson loop operator. Repeating the steps from 
Eq.qnjD to Eq.(]T3[) we get for < W{C t ) > 

<W(C t )>=^(2nP v )- N ^ £ exp (__L< n , n> ). (19) 

n 

(dn = C t ) 

This is a sum over all world-sheets spanned on the loop Ct- In other words, we have arrived 
to an expression of the partition function of compact electrodynamics in terms of the world- 
sheets of pure electric loop excitations. Notice that the loop action corresponding to the 
Villain action is proportional to the quadratic area A 2 : 

i.e. the sum of the squares of the multiplicities s p of plaquettes which constitute the loop's 
world sheet S. It is interesting to note the similarity of this action with the continuous 
Nambu action or its lattice version, the Weingarten action p3 , which is proportional to the 
area swept out by the bosonic string. 



III. NUMERICAL ANALYSIS 



We simulated the loop lattice action of eq. (pup. In ref. (TB| we shown that this action is 
equivalent to the Integer Gaussian Model (IGM) that, in turn, it can be also reached as a 
limit of the non compact scalar QED in the limit of infinite Higgs coupling. In that paper we 
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presented a direct simulation of the IGM which from the analysis of the energy histograms 
showed clearly the first order nature of the transition. Nevertheless, we did not computed 
the critical exponents of the phase transition. 

Here we have performed a more quantitative study and computed the exponents of the 
phase transition. We applied the simulation method introduced in ref. |p.7[ . This method 
consists in generating collections of integer variables attached to plaquettes which form 
surfaces with the appropriate constraints for each model (null frontier in the present pure 
gauge case). The variables are updated by means of a standard Metropolis Monte Carlo 
algorithm. These larger simulations produced a phase transition for the coupling constant 
(3 = l/(2/V) = 0.778 i.e. (3 V = 1/(2 * 0.778) = 0.6426 which is slightly higher than the 
previous value of j3y — 0.64 |p7| . 

The physical magnitudes measured under the simulation are the internal energy and the 
specific heat both normalized to the number of plaquettes N p = 6L 4 : 

e = N- 1 (E) (21) 
C v = (3'N- 1 ((E 2 ) - (E) 2 ) (22) 

on (3 = 1/2/V- 

The details of the run are the following: 

• Hypercubical lattice (that is, periodic boundary conditions) 

• Sizes ranging from L = 8 to L = 24. 

• Good statistics, more than 100 times the self-correlations. It means millions of itera- 
tions for L = 24. 

• In order to locate the transition and the specific heat peaks, we didn't do thermal 
cycles. Instead we used iteratively Ferrenberg-Swedsen reweighting until a two peaks 
histogram were found. 

We made a FSS analysis for the following quantities: 

<0 The specific heat maxima C max (L) (§ [IV A ). 



<0 The critical coupling position f3 c (L) (§ |1V B| ). 

(} The partition function zeroes in the complex plane. (§ |IV C|) . 

<C> The latent heat, that is the distance between the two peaks (§[L V Dp. 



The calculations have been performed in Pentium computers running Linux (including 
the RTNN machine of 32 Pentium Pro processors at 200MHz). As it is expected, the elapsed 
time in the runs grows with the volume of the lattice obeying the following relation: 

{ L \ 4 

t = 42 I — J hours/Mit, (23) 



Vio 

where t is the elapsed time for a Pentium Pro 200MHz running a million of iterations 
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We have been very careful with the statistical analysis of such a transition. We have 
measured the autocorrelations, r e , of the time-series for each lattice size in order to determine 
both, the number of thermalization iterations and the number of measures to obtain a 
significative number of independent measures. It is accepted [|4| that about 20r e iterations 
are enough in order to have the series thermalized and about 100r e to have a sufficiently large 
time series. Table | collects this data. Using this data, one obtains that the self-correlation 



L 


Pmc 


r 


nth/r 


n m /r 


6 


0.77991 


90 


1100 


110000 


8 


0.77821 


300 


330 


660 


10 


0.77821 


500 


100 


200 


12 


0.77700 


810 


70 


140 


14 


0.77700 


1500 


60 


100 


16 


0.77688 


2100 


50 


225 


18 


0.77688 


3300 


40 


130 


20 


0.77688 


5000 


30 


110 


24 


0.77680 


13100 


20 


100 


TABLE I. Correlations 



time grows exponentially (critical slowing down) according to 

T = T e 2aL r = 49(4) a = 0.116(2) Q = 0.95, (24) 

where we supposed an error of about a 5% measuring r.Q 
The goodness of the fits is done using 

where Q(a,x) is the regularized gamma function 

T(a,x) 1 r°° 



Q( a? x ) = i- P( a , x) = ee — - I e-H^dt. (26) 

1 (a) 1 {a) Jx 



Combining (P3|) and ( f24|) we obtain the CPU times needed for an arbitrary simulation. 
We spent about three months of calculations for L = 24, and a simple estimation reveals 
that to calculate L = 28 in our machine we would need more than a year. 

In Fig. [I] one can see energy histograms for different lattice sizes an the time series 
showing the presence of tunneling between the two phases. This undoubtedly shows that 
the world-sheet formulation (in terms of integer variables on plaquettes) is considerable 
more effective than the standard one (in terms of gauge potentials on links) to carry out a 
simulation. 



It is difficult to give a realistic estimation for the error in r. Nevertheless, even underestimating 
the errors (1%) the goodness of the fit is not bad at all (Q = 0.04). 
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0.26 - 




200 400 0.03 0.06 

10 iteracions P(e) 

FIG. 1. Some Monte Carlo runs with their corresponding histograms for L = 12,16,20,24. 
All the axes are equally scaled in order to compare them. Notice that for large L the separation 
between phases grows, that is, the two peaks structure is more clear and in the time-series, the 
chance for tunneling decreases. 
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IV. CRITICAL EXPONENTS 

Table |I| summarizes all the observables measured in this work. 



L 


Pmc 






lO- 3 • X (°) 


°« 


Ae 


6 


0.77991 


0.7806(5) 


0.780215(2) 


9.9(2) 


3.4750(7) 




8 


0.77821 


0.7785(4) 


0.778358(3) 


4.3(1.8) 


5.42(6) 




10 


0.77821 


0.7780(2) 


0.777978(4) 


2.0(2) 


8.44(9) 


0.0237(5) 


12 


0.77700 


0.7770(2) 


0.776985(1) 


1.2(2) 


10.2(3) 


0.0206(4) 


14 


0.77700 


0.7769(7) 


0.776972(8) 


0.76(2) 


14.7(4) 


0.019(1) 


16 


0.77699 


0.7769(5) 






20.1(5) 


0.0169(1) 


18 


0.77688 


0.7768(2) 


0.776843(2) 


0.36(3) 


24.2(2) 


0.0155(2) 


20 


0.77688 


0.7768(8) 






36.2(7) 


0.0160(5) 


24 


0.77680 


0.7767(7) 




0.23(2) 


65.8(1.8) 


0.0145(4) 



TABLE II. Observables measured in this work. 1st an 2nd column contain the size of the 
lattice and Pmc of the run. 3rd column contains the value of the coupling, C , of the maximum of 
C v . In 4th and 5th columns one can find the real and imaginary part of the coupling that vanishes 
Z. 6th column contains C v and finally 7th column is the latent heat, Ae. 



A. Specific heat 

In a 1st order phase transition, the specific heat is expected to have a volume scaling 
law, i.e. 

C™ ax (L) = a l + b l L D . (27) 

But a continuous phase transition must obey 

C™ ax {L) = a 2 + b 2 L a ' v . (28) 

Both expressions are equal if a = 1 and v — 1/D. Using the data collected in column 6 of 
the table || we have performed a fit according to those scaling laws. Results are summarized 
in table 

Fig. shows the measured values of C™ ax as a function of L 4 . Notice the alignment of 
the points, exhibiting a 1st order behavior. Solid line is the 1st order fit using only values 
14 < L < 24. 

B. Critical coupling 

The critical coupling must reach their asymptotic behavior obeying the following laws, 
P e (L) = &(oo) + AL~ 1/U (2nd order), (3 C (L) = (3 c (oo) + AL~ D (1st order). (29) 
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1st order, eq. (|27|) 








2nd order, eq. 


(H) 


L 


a\ 


10" 4 6i 


X 2 /d.o.f 




Q 




a/u 


X 2 / d.o.f 


Q 


6-24 


3.16(6) 


2.5(3) 


2.9 • 10 2 


0.51 


■ 10~ 


-375 


2.44(4) 


18.8 


0.97 • 10 


8-24 


5.08(5) 


2.1(3) 


7.3 • 10 


0.38 


• 10" 


-76 


2.7(7) 


19.2 


0.91 • 10~ 15 


10-24 


6.7(1) 


1.9(3) 


5.1 


0.43 


• 10" 


-3 


3.60(13) 


1.5 


0.22 


12-24 


7.7(3) 


1.8(3) 


2.8 


0.38 


• io- 


-1 


3.4(2) 


1.6 


0.19 


14-24 


8.6(4) 


1.7(5) 


5.8 • 10~ 2 


0.94 






3.9(3) 


0.079 


0.78 


16-24 


8.8(9) 


1.7(6) 


6.2 • 10" 3 


0.94 













TABLE III. Fits of C™ ax . The 1st order fit was done using eq. (|27|) (growing proportional to 
L D ). 2nd order fit corresponds to eq. d28|), that is, growing as L a l u , being a/v free. Notice that 
in both cases the best fits are obtained using only the values 14 < L < 24. In addition, the best 
2nd order fit gives and exponent close to 4. 



It is common to say that 1st order transition have a "critical exponent" of v = 1/D (=0.25 
in 4 dimensions). Trivial 2nd order transitions have v = 0.5. Nevertheless, some weak 1st 
order transitions have an intermediate exponent and is not easy to elucidate the order of 
the transition from this scaling law. In our case, the errors are so large that all the fits are 
equally "good" and therefore this quantity does not provide much evidence on the order of 
the phase transition. 

C. Lee- Yang zeroes 

The zeroes of the partition function analytically extended to the complex plane manifest 
a scaling behavior governed by the critical exponent — 1/v p5|-|28[] Table show the location 



of the zeroes closest to the real axis. From this data, a fit to the imaginary part of (3 gives 

X®(L) = BL-i (30) 

This quantity is rather insensitive to the finite size effects. The best fit (5 < L < 20) 
gives 

B = 2.3(2) v = 0.329(3) (31) 

with x 2 / d.o.f = 0.6 and Q = 0.7. Fig |3] show the data points and the regression curve. 

D. Latent heat 

The last column of table [TI| show the latent heat for each lattice size. We expect a 
behavior as 

Ae(L) = Aei(oo) + ciiT* (32) 



10 



80 




N/6 

FIG. 2. FSS of the specific heat as a function of the volume. The points with error bars are 
the data and the solid line is the fit of eq. ( |27| ) using the three last points. Notice that the smaller 
lattices are away from this fit. 

This gives the result Aei(oo) = 0.0117(13) with d x = 1.8(4). {x 2 /d.o.f = 0.8 i Q = 0.53) 
We have also performed a fit as 

Ae(L) = Ae 2 (oo) + c 2 e- d2L (33) 

with the result Ae 2 (oo) = 0.0135(6) i d 2 = 0.19(3) (x 2 /d.o.f = 0.5 i Q = 0.75) 

Figure ^ show these two different fits together with the measured data. In both cases 
and asymptotic behavior with Ae = is ruled out. 

V. CONCLUSIONS 

Among different appealing features the world-sheet formulation avoids the summation 
over redundant gauge configurations. This fact implies an amelioration of the convergence 
to the equilibrium. As a consequence, the time series generated using this approach show a 
high degree of "tunneling" between phases, a phenomenon that is depressed in the standard 
Wilson formulation. 

We have observed in the simulations a standard first order behavior in several observables 
(autocorrelations, specific heat, latent heat) with an intermediate behavior between the 
first order [y = 0.25) and second order Gaussian exponent [y = 0.5) like other weak first 
order transitions. Wether or not this observed first order signature is an artifact of the 
topology/boundary-conditions of the lattice is an issue we will analyze in a future work. 
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10 4 I 1 1 1 1 1 1 1 1 

6 8 10 12 14 16 18 20 

L 



FIG. 3. FSS of x^°\ the imaginary part of the zeroes in the partition function. This is a log-log 
plot, so the fit must be a line (the solid line) 



0.025 









\ 


-d L 

Ae(L) = Ae 2 +c 2 e 2 




\ 


Ae(L) = Ae l +c 1 L~ c ' 1 






- Ae, = 0.01 17(13) 






Ae 2 = 0.0135(6) 
















I.I.I 



10 15 20 25 30 

L 

FIG. 4. FSS of the latent heat. Data fit so well an exponential law p^ ) and a potential one (|32|). 
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